TWO-BODY DIRAC EQUATIONS FOR 
RELATIVISTIC BOUND STATES OF QUANTUM FIELD THEORY 

Horace W. Crater 
The University of Tennessee Space Institute 
Tullahoma, Tennessee 37388 
Peter Van Alstine 
12474 Sunny Glen Drive, Moorpark, Ca. 93021 

Abstract 



We review a little-known treatment of the relativistic two-body bound- 
state problem - that provided by Two-Body Dirac Equations obtained from 
constraint dynamics. We describe some of its more important results, its 
relation to older formulations and to quantum field theory. We list a number 
of features crucial for the success of such a formulation, many of which are 
missing from other methods; we show how the treatment provided by Two- 
Body Dirac Equations encompasses each of them. 



Invited paper presented at a Conference on September 12th, 1997 at the University of 
Georgia in honor of Professor Donald Robson on his 60th birthday. 
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I. INTRODUCTION 



The correct formulation of relativistic two-body bound state wave equations and their 
connection to quantum field theory is an old problem going back to papers by Eddington and 
Gaunt |]J in 1928 (bases on conjectures by Heisenberg). But judging from the large variety 
of approaches attempted even in recent years, this problem has no generally agreed-upon 
solution. Perhaps for this reason, recent field theory books have avoided this topic. For 
example, Steven Weinberg || states: "It must be said that the theory of relativistic effects 
and radiative corrections in bound states is not yet in satisfactory shape." 

In H we have listed many of the attempts at solving this problem made over the past 
47 years. At the head of this list, appear the work of Nambu and the better- known work of 
Bethe and Salpeter, both developed over 20 years after the partly successful, semirelativistic 
equation of Breit. The Bethe-Salpeter equation is an integral equation in momentum space 
that is manifestly covariant, obtained directly from relativistic quantum field theory. Over 
the years, however, many problems have turned up to impede its direct implementation. 
These are the sources of the numerous attempts to reformulate the two body problem of 
relativistic quantum field theory. For example as pointed out by Weinberg || : "The 
uncrossed ladders can be summed by solving an integral equation known as the Bethe- 
Salpeter equation, but there is no rationale for selecting out this subset of diagrams unless 
both particles are non-relativistic, in which case the Bethe-Salpeter equation reduces to 
the ordinary nonrelativistic Schrodinger equation, plus relativistic corrections associated 
with the spin-orbit couplings that can be treated as small perturbations." Furthermore, the 
Bethe-Salpeter equation in the ladder approximation possesses negative norm or ghost states, 
due to its treatment of the relative time degree of freedom - spoiling the naive interpretation 
of it as a quantum wave equation. 

Salpeter and many others p| have developed noncovariant instantaneous truncations of 
this equation. For a variety of reasons most of the attempts we have cited are not appropriate 
for the treatment of highly relativistic effects like those necessary for the calculation of quark- 
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antiquark bound states. 

Among the authors who have tried to rectify this problem are Professor Robson and 
his collaborator Dr. Stanley who gave one of the first comprehensive attempts to obtain 
a potential-model description of the entire meson spectrum, combining exact relativistic 
kinematics with a non-perturbative treatment of the effects of the important spin- dependent 
interactions. On the other hand, our work treats much the same spectral phenomenology 
but remains close to Dirac's own one body work by starting with a pair of covariant two- 
body Dirac equations (one for each particle) which forces certain restrictions on the various 
spin-dependent interactions that can appear. 

II. TWO-BODY DIRAC EQUATIONS FROM CONSTRAINT DYNAMICS 

The list of references on the relativistic two-body problem includes diverse alternative 
approaches. On the other hand the one-body Dirac equation Q has no serious rivals. That 
equation is a well defined wave equation that can be solved nonperturbatively and serves as 
an example of a successful bound state equation. 



provides a relativistic version of Newton's first law, with associated relativistic second law 
appearing as the four-vector substitution 



A. The One-Body Dirac Equation 



The free Dirac equation 



(7 • p + m)ip = 



(2.1) 



(2.2) 



for electromagnetic interaction and as the minimal mass substitution 



m 



m + S 



(2.3) 
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for scalar interactions combining to give 

( 7 - (p- A) +m + S)i) = 0.2.4 (2.4) 

As we shall see, the two-body Dirac equations take simple forms that generalize this one to 
the interacting two-body system. 



B. The Breit and Eddington-Gaunt Equations 

The Breit equation || and most three dimensional truncations of the Bethe-Salpeter 
equation are neither well defined wave equations, nor manifestly covariant. Breit proposed 
his equation equation 



Em 



|«i • pi + ftrnj + a 2 ■ p 2 + fornv - "[1 - -(ai • a 2 + ai • r a 2 • r)]j \P (2.5) 



in 1929 as a correction to an earlier defective equation which he called the Eddington-Gaunt 
equation. 

The Eddington-Gaunt equation |1[ has the form 

E^) = jai • pi + (3 1 m 1 + a 2 -p 2 + - ~(1 - «i • "2)| (2.6) 

This equation was defective because it failed to include the semirelativistic electrodynamic 
interaction of Darwin ||, together with the Coulomb interaction in the combination 

a 1 

1 - -(vt ■ v 2 + V! ■ rv 2 ■ f) . (2.7) 

r 2 

This structure arises from the semirelativistic expansion of the Green function, including 
retardative terms through order 4-. 



J J J\GJ 2 = J dr\ J dT 2 XiX2 fM 5[(xi - x 2 ) 2 } 



a 1 

— 1 - ~(vi ■ v 2 + vi ■ f v 2 ■ f)\. (2.8) 
r 2 
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Replacing velocities in the Darwin interaction by dps led Breit to terms missing from Ed- 
dington's approach. 

In modern language the Eddington-Gaunt interaction is most closely connected to QED 
in the Feynman gauge while the Breit interaction is most closely connected to QED in the 
Coulomb gauge. Eddington's incorrect implementation of the Feynman gauge produces the 
wrong QED spectrum while Breit's correct implementation of the Coulomb gauge produces 
the correct QED spectrum for positronium and muonium. 

Since the Breit equation is not covariant and not a well defined wave equation it must 
be handled by semirelativistic perturbative methods. For example, for QED it is standardly 
rearranged as the semirelativistic elaboration of the Schrodinger Equation: 

Hip = wip. (2.9) 

(w is the total cm. energy) in which the Hamiltonian is the first-order perturbative form 

2m x 8mf 2m 2 8m 2 r 

p 2 1 1 

-<*( - [ h V ■ (1 - ft) ■ p\ordered 

m\mi r 2m 2 m 2 r 



1,1 1 1 L 1 2 _ .1 2 w . 

- — + —)6(r) - -- ■ — + Vi + — + — y 2 

2 mf wij 4 r 6 mf mim 2 m\m2 



-a 1 -a 2 5{r) + = )). (2.10) 



Am\m 2 3 r 3 r 5 



C. Manifestly Covariant Two-Body Dirac Equations 

In the 1970's, several authors used Dirac's constraint mechanics to attack the relativistic 
two-body problem at its classical roots 0] successfully evading the so-called no interaction 
theorem 0. Using this method, the present authors extended the constraint approach to 
pairs of spin one half particles to obtain two-body quantum bound state equations that 
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correct || the defects in the Breit equation and more importantly in the ladder approxima- 
tion to the Bethe-Salpeter equation, exorcising quantum ghosts by covariantly controlling 
the relative time variable. Those equations are the two-body Dirac equations of constraint 



which correct defects in patchwork approaches. They are manifestly covariant while yield- 
ing simple three-dimensional Schrodinger-like forms similar to those of their nonrelativistic 
counterparts. Their spin dependence is not put in by hand, as in patchwork approaches, but 
is determined naturally by the Dirac-like structure of the equations. These equations have 
passed numerous tests showing that they reproduce correct QED perturbative results when 
solved nonperturbatively. They thus qualify as bona fide wave equations. In addition, the 
Dirac forms of the equations automatically make unnecessary the ad hoc introduction of cut- 
off parameters. The relativistic potentials appearing in these equations are related directly 
to the interactions of perturbative quantum field theory or (for QCD) may be introduced 
semiphenomenologically. In QCD one can regard these manifestly covariant equations as an 
anticipation of those that may eventually emerge from lattice gauge theory as applied to 
meson spectroscopy. 

As in the ordinary one-body Dirac equation Eq. Q2.4p , for particles interacting through 
world vector and scalar interactions the two-body Dirac equations equations take the general 
minimal coupling form 



These equations provide a non-perturbative framework for extrapolating perturbative field 
theoretic results into the highly relativistic regime of bound light quarks, in a quantum 
mechanically well defined way. That framework incorporates of three main properties: 
I: Exact Lorentz covariance, 
II: Minimal interaction structure 



dynamics ||10|| . They possess a number of important features some of which are unique and 




(2.11a) 




(2.11b) 
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Ill: Compatibility of the two-equations (which leads to a relativistic 3rd law, covariantly 
restricting the relative momentum and energy while correctly structuring spin-dependent 
recoil) : 



[5i,5 2 ]^ = 0. 



(2.12) 



The satisfaction of this requirement originates in part from the presence of supersymmetries 
in the (pseudo-) classical limit of the Two-Body Dirac Equations. The potentials in these 
equations for four-vector and world scalar interactions are intimately connected to those of 
Wheeler- Feynman Electrodynamics (and its scalar counterpart) |TT|JT0| , p!2[| and have been 
obtained systematically from perturbative quantum field theory fll3"| , |14j| . 

For vector interactions alone, their momentum and spin dependences take the simple 



"hyperbolic" forms p5|JT6[] 



Ai = [1 - cosh(£)]pi + sinh(£)p 2 - ~{de G ■ 72)72 



(2.13a) 



in which 



A 2 = [1 - cosh(£)]p 2 + sinh(£)pi + -(de s ■ 71)71 



(2.13b) 



g = -^H 1 - ia/w) 



(2.14) 



(with w the total cm. energy). We originally found the logarithm form in a derivation from 
the Wheeler- Feynman classical field theory. In fact, in quantum electrodynamics that form 
turns out to embody an eikonal summation of ladder and cross-ladder diagrams. 

In Eqs. (|2.14j ), the invariant A is a function of the covariant spacelike particle separation 



= x » + p^(p . x ) 



(2.15) 



perpendicular to the total four-momentum, P. (P = £ is a time-like unit vector.) Its 
appearance signifies that the dynamics is independent of the relative time in the cm system. 
For lowest order electrodynamics, 



A=A(x±) = -- (2.16) 
r 

in which 

r = \fx\. (2.17) 

The form of the covariant spin-dependent terms and the fact that A depends on x± are 
consequences of compatibility of the two Dirac equations ([Si,^]^ = 0). In quark- model 
calculations, the invariant A and its counterpart for the scalar interaction are chosen on 
semiphenomenological grounds. 

These two-body Dirac equations bypass most of the difficulties of the Bethe-Salpeter 
equation that arise from the presence relative time and energy variables, and yield a three- 
dimensional but manifestly covariant rearrangement of the Bethe-Salpeter equation. The 
three dimensional character is partially embodied in the invariant r which reduces to the 
interparticle separation only in the cm. system. The fact that the interaction is instanta- 
neous in the cm. system is a direct consequence of the compatibility of the two equations. 
It is not an ad hoc restriction imposed on the equation as is done for various instantaneous 



approximations of the Bethe-Salpeter equation [|1J . 

Just as Dirac arrived at his equation by "taking the square root of the Klein-Gordon 
equation" so these equations can be derived by rigorously "taking the square root" of the 
corresponding compatible "two body Klein-Gordon equations" [O 



S 1 = 751(71 • (Pi - A x ) + mi + = "J (Pi - A) 2 + K + Sx) 2 + ..." (2.18a) 



S 2 = 7si(7i ■ (Pi - Ax) +mx + Sx) = V (Pi - Ax) 2 + (mx + Sx) 2 + ...". (2.18b) 



It is in this sense that it is most natural to call Eqs.( |2.1lD "Two-Body Dirac equations 



These equations are not only manifestly covariant but are also quantum-mechanically 
well-defined. That is, their covariant Schrodinger-like forms for the effective particle of 
relative motion 
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(2.19) 



in which 



6 2 H 



: (w — 2(m 1 + rro 2 )io + (m 1 



2 2\2 



2) )i 



(2.20) 



can be solved nonperturbatively for both QED and QCD bound state calculations since 
every term in $ w ((Ti, 02, -A(r), S(r)) is well defined (less singular than — l/4r 2 ) |T7|JT4]]. Fur- 



thermore, recent work has shown that the Schrodinger-like forms can be transformed into 
equations that, like their nonrelativistic counterparts, involve at most 2 coupled wave func- 
tions |I8|JT9[| even when non-central tensor forces or spin-difference-orbit interactions are 
present. Note that in the decoupled form (or any other convenient form) the specific forms 
of the spin dependent potentials are dictated (through the reduction process) by the inter- 
action structure of the original Dirac equation Eq. [2.11| and are not put in by hand. 

We have checked the nonperturbative validity of these equations as well defined wave 
equations by solving them analytically and numerically to obtain the standard fine and 
hyperfine spectra of QED. For example, we obtained an exact spectral solution for the 
singlet positronium system A = —a/r 



w = m 



i+ 



TO 



1)2 



a' 



(2.21) 



of the fully coupled system of 16-component equations 15,20 



Slip = S21/) = 0. 



Such validation we claim ought to be required of all candidate equations for nonperturbative 
quark model calculations. No others on that list of 60, including the Salpeter reduction (the 
no-pair Breit equation) and the Blankenbecler-Sugar equation, have yet met this demand. 
If not required to meet this demand two body formalisms may lead to possibly spurious 
nonperturbative predictions. For example Spence and Vary ]22| have shown that in their 
treatment, the no-pair Breit equation and the Blankenbecler-Sugar equation predict low 
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energy electron-positron resonances between 1.4 and 2.2 MeV in Bhabba scattering. Such 
states have not been observed in low energy electron-positron collisions nor are they predicted 
by our equations. With no rigorous check on nonperturbative solutions for ordinary QED 
bound states, how can the QCD results of such treatments be trusted? 

Given a static potential model V = V(r) for the quark- ant iquark interaction we can 
incorporate it in a covariant way into our equations by ||i~7| . |T4| 

a) replacing nonrelativistic r by 

b) parcelling out the static potential V into the invariant functions A(r) and S(r). This 
step remains a partially phenomenological one. However, in our approach once A and S are 
fixed so are all the accompanying spin dependences. One cannot adjust the various parts 
independently, as is done in many approaches [^1[ which just add plausible potential energy 
terms to two-body relativistic kinetic energy operators. 



D. Comparison of Two-Body Dirac with Breit, Eddington- Gaunt and other 

Approaches 

Since the Two-Body Dirac Equations make quantum mechanical sense they can be solved 
nonperturbatively in QCD bound state calculations, just like the ordinary one-body Dirac 
equation (to which they reduce in the limit that either particle becomes infinitely massive. 
The well-defined potential structures function as a natural smoothing mechanism that yields 
the correct spectrum while avoiding singular effective potentials like delta functions that 
appear in the Breit reductions and most competing approaches. For example, the Pauli, 
i.e. Schrodinger-like forms of our equations make quantum mechanical sense in the strong 
potential, nonperturbative regime where relativistic effects of the wave operator on the wave 



function are not negligible. This claim is easiest to see by examining the connection |T7| of 
the main spin-spin term in our equation with that of Breit. 

(Two - Body Dirac) - ■ a 2 d 2 ln(l - — ) -> - ai ' °2^A Breit (2.22) 

6 W 3 mi + 777.2 
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For ,4's that have singular short range behaviors like —a/r (QED) and 8ir/27rlnr (QCD) 
the weak potential form which appears in the reduced Breit equation and most patchwork 
approaches can only be used in a perturbative calculation. The two-body Dirac form on the 
left can be used nonperturbatively with the logarithm term providing a natural smoothing 
mechanism avoiding the necessity of introducing singularity softening parameters in phe- 
nomeno logical approaches [jl7| ,[l4}|. 



For QED, the Two-Body Dirac equations work naturally in the covariant Feynman gauge, 
organizing the diagrammatic summation in a more efficient way than does the BSE. They 
naturally produce an equation closely connected in form to the defective equation of Ed- 
dington but yielding a correct spectrum. They achieve this result by effectively summing an 



infinite number of ladder and cross ladder diagrams in a kind of eikonal approximation p3 
Can we see this effect directly? Since they form a compatible pair they can be combined in 
any number of equivalent ways, in particular in a Breit-like form |T^j9|,p^.p4 



w \fr = {a x ■ px + f3imi + a 2 ■ p 2 + l3 2 m 2 + w(l - exp [ - G(x ± )(l - 3.x • a 2 )])} V , (2.23) 

1 2r> n 

g= --ln(l + —) = -— + ... (2.24) 

2 rw wr 

Comparing this to Eddington and Gaunt 's 

= j<3i ■ Pi + Ami + 6l 2 -p 2 + (3 2 m 2 - "(1 - ol x ■ 5 2 ) X I / | , (2.25) 

we see that by effectively stopping at lowest order, the Eddington-Gaunt equation was 
doomed to failure. 

In detail, to all orders in the potential, our equations produce 



w{l - exp[-g(x±)(l -ai -a 2 )]) (2.26) 

A 2 A 3 1 

= ^.(1 - ax ■ a 2 ) (1 - 3x ■ o 2 ) (1 - 751752 + «i • a 2 - ox ■ a 2 ) . (2.27) 

w w z {1 — zA/w) 

^From the point of view of the Breit form of our equations this means that not only do they 

contain vector interactions but other covariant interactions as well (e.g. pseudovector, scalar 
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and pseudoscalar interactions), which plausibly can be viewed as originating in the sorts of 
products of vector interactions that occur in multiparticle exchange. Using a transformation 
due to Schwinger, we have shown elsewhere that the extra terms that remain in the weak 



potential limit the are canonically equivalent to Breit's retardative terms [16 



E. Connection to Quantum Field Theory 



The invariant forms A and S in the Two-Body Dirac Equations Eq. fl2.llp may be sys- 



tematically obtained from the corresponding quantum field theories |]13|,^4| . The connection 
is 

(7 2 , A(r), S(r)) = m5(P ■ p)K{l + IC)' 1 (2.28) 
giving the quasipotential in terms of the Bethe-Salpeter kernel /C and its projection 

£ = GJC (2.29) 

in which 

1 1 

G = (-= = --= = - m5(P ■ p)- 2 — -) (2.30) 

K pj + ml - tO pi + ml - tO v ' p\ - b 2 (w) - t0 J v ' 

is the difference between forms of two-body propagators as given by the Bethe-Salpeter 
equation on the one hand and the constraint equations on the other (with the relativistic 
third law delta function.) (p is the relative momentum and P is the total momentum.) 
This connection is a sort of covariant version of the three-dimensional Lippman-Schwinger 
equation. The difference is that in this equation, derived by Sazdjian |L3| as "quantum 
mechanical transform of the Bethe-Salpeter equation" , the potential follows from the irre- 
ducible scattering matrix rather than the other way around as done in the nonrelativistic 
case. The present authors have also derived the effective potentials from classical field the- 
ory in the form of Fokker- Tetrode actions through comparison with the classical limit of the 
constraint equations [P^,P~3 • 
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III. TWO-BODY DIRAC EQUATIONS AND MESON SPECTROSCOPY 



We have formulated a constraint version of the naive quark model for mesons by using the 
static Adler-Piran quark-antiquark potential [|26|| , (covariantly reinterpreted) in our Two- 
Body Dirac Equations. Adler and Piran obtained their static quark potential from an 
effective non-linear field theory derived from QCD. It has the general form 

V AP (r) = A(U(Ar) + U ) (=A + S). (3.1) 

Since their potential is nonrelativistic it cannot distinguish between world scalar and vec- 
tor potentials, simply representing the effect of their sum in the nonrelativistic limit. It 
incorporates asymptotic freedom analytically through 

AC/(Ar « 1) ~ l/(rlnAr) (3.2) 

and linear confinement through 

AU(Ar » 1) ~ A 2 r. (3.3) 

At long distances their potential includes not only the linear confinement piece but also 
subdominant logarithm terms among others 

V A p(r) = A(dx + c 2 ln(x) + c 3 /Vz + c 4 /x + c 5 ), x = Ar > 2. (3.4) 

The CiS are given by the Adler-Piran leading log-log model. The realistic Adler-Piran po- 



tential or ones like it such as the Richardson potential [27], fail miserably for light mesons 



when used in the nonrelativistic Schrodinger equation [28]. Exact covariance is essential to 
handle the light mesons if one insists, as we do, on using potentials closely tied to QCD. 

We apportion the potential between the relativistic invariants S and A that determine 
the scalar and vector potentials according to the scheme 

A = exp{-/3r)[V AP - -] + - + (3.5a) 

rr> rr> rr> 

S = V AP + — -A. (3.5b) 
r 
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In this way we impose the requirement that at short distance the potential is strictly vector 
while at long distance the confining portion is scalar with Coulomb vector portions. The 
relativistic invariance of S and A follows by reinterpreting the variable r as r = 

Our quark model is a naive quark model in that we ignore flavor mixing and the effects 
on the bound state energies of decays. The results we obtain, using the same potentials for 
all of the mesons, are spectrally quite accurate, from the heaviest upsilonium states to the 
lowly pion. 

TABLE I - MESON MASSES FROM CO VARIANT CONSTRAINT TWO- 
BODY DIRAC EQUATIONS 

NAME EXP. THEORY 



T : 


bb 1 3 S 1 


9.460( 


0.2) 


9.453( 


0.6 


T : 


bb 1 3 P 


9.860( 


1.3) 


9.842( 


1.4 


T : 


bb 1 3 P 1 


9.892( 


0.7) 


9.889( 


0.1 


T : 


bb 1 3 P 2 


9.913( 


0.6) 


9.921( 


0.5 


T : 


bb 2 3 Si 


10.023( 


0.3) 


10.022( 


0.0 


T : 


bb 2 3 P 


10.232( 


0.6) 


10.227( 


0.2 


T : 


bb 2 3 P l 


10.255( 


0.5) 


10.257( 


0.0 


T : 


bb 2 3 P 2 


10.269( 


0.4) 


10.277( 


0.8 


T : 


bb 3 3 Si 


10.355( 


0.5) 


10.359( 


0.1 


T : 


bb 4 3 Si 


10.580( 


3.5) 


10.614( 


0.9 


T : 


bb 5 3 S 1 


10.865( 


8.0) 


10.826( 


0.2 


T : 


bb 6 3 S l 


11.019( 


8.0) 


11.013( 


0.0 


B : 


bu l 1 ^ 


5.279( 


1.8) 


5.273( 


0.1 


B : 


bd l 1 ^ 


5.279( 


1.8) 


5.274( 


0.1 


B* 


: bu l 3 ^! 


5.325( 


1.8) 


5.321( 


0.1 


B s 


: bs I 1 Sq 


5.369( 


2.0) 


5.368( 


0.0 


B s 


: 6s l 3 ^ 


5.416( 


3.3) 


5.427( 


0.1 
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Vc ■ 


cc 1% 


2.980( 


2.1) 


2.978( 0.0) 


V: 


CC l 3 Si 


3.097( 


0.0) 


3.129( 12.6) 


Xo 


CC l X Pi 


3.526( 


0.2) 


3.520( 


0.4) 


Xo 


cc 1 3 P 


3.415( 


1.0) 


3.407( 


0.4) 


Xi 


CC l 3 Pi 


3.510( 


0.1) 


3.507( 


0.2) 


X2 


CC 1 3 P 2 


3.556( 


0.1) 


3.549( 


0.6) 


Vc ■ 


cc 2% 


3.594( 


5.0) 


3.610( 


0.1) 


^: 


cc 2 3 Si 


3.686( 


0.1) 


3.688( 


0.1) 


V>: 


CC 1 3 P>! 


3.770( 


2.5) 


3.808( 


2.0) 


^: 


cc 3 3 Si 


4.040( 10.0) 


4.08K 


0.2) 


^: 


CC 2 3 P! 


4.159( 20.0) 


4.157( 


0.0) 


^: 


CC 3 3 P>! 


4.415( 


6.0) 


4.454( 


0.4) 


D : 


cu l 1 ^ 


1.865( 


0.5) 


1.866( 


0.0) 


£> : 


cd l 1 ^ 


1.869( 


0.5) 


1.873( 


0.1) 


P>* 


: cu l 3 ^! 


2.007( 


0.5) 


2.000( 


0.4) 


P>* 


: erf l 3 Si 


2.010( 


0.5) 


2.005( 


0.3) 


D* 


: ell l 3 Pi 


2.422( 


1.8) 


2.407( 


0.6) 


D* 


: cd l 3 Pi 


2.428( 


1.8) 


2.411( 


0.5) 


D* 


: ell 1 3 P 2 


2.459( 


2.0) 


2.382( 


11.3) 


D* 


: cd 1 3 P 2 


2.459( 


4.0) 


2.386( 


3.5) 


D s 


: cs l 1 ^ 


1.968( 


0.6) 


1.976( 


0.5) 


D* s 


: cs 1 3 S 1 


2.112( 


0.7) 


2.123( 


0.9) 


D* s 


: cs l 3 Pi 


2.535( 


0.3) 


2.511( 


6.2) 


D* s 


: cs 1 3 P 2 


2.574( 


1.7) 


2.514( 


9.6) 


K 


su l 1 ^ 


0.494( 


0.0) 


0.492( 


0.0) 


K 


sd l 1 ^ 


0.498( 


0.0) 


0.492( 


0.4) 


K* 


: su l 3 Si 


0.892( 


0.2) 


0.910( 


0.6) 


K* 


: sd l 3 ^ 


0.896( 


0.3) 


0.910( 


0.3) 
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: su 


1.273( 7.0) 


1.408( 


3.2 


K* 


: sU 1 3 P 


1.429( 4.0) 


1.314( 


0.7 


K x 


: l 3 Pi 


1.402( 7.0) 


1.506( 


1.0 


K* 


: sU 1 3 P 2 


1.425( 1.3) 


1.394( 


0.5 


K* 


: srf 1 3 P 2 


1.432( 1.3) 


1.394( 


0.6 


K* 


: su 2 1 So 


1.460( 30.0) 


1.591( 


0.2 


K* 


: su 2 3 S! 


1.412( 12.0) 


1.800( 


6.7 


K 2 


: 1 X P 2 


1.773( 8.0) 


1.877( 


0.8 


K* 


: su l 3 P>i 


1.714( 20.0) 


1.985( 


1.4 


K 2 


: SM 1 3 P> 2 


1.816( 10.0) 


1.945( 


1.3 


K 3 


: su 1 3 P 3 


1.770( 10.0) 


1.768( 


0.0 


K* 


: su 3 1 S 


1.830( 30.0) 


2.183( 


1.4 


K* 


: sU 2 3 P 2 


1.975( 22.0) 


2.098( 


0.2 


K* 4 


: su 1 3 P 4 


2.045( 9.0) 


2.078( 


0.1 


K 2 


: sm 2 3 P> 2 


2.247( 17.0) 


2.373( 


0.5 


K* 


: su 1 3 G 5 


2.382( 33.0) 


2.344( 


0.0 


K* 


: su 2 3 P 3 


2.324( 24.0) 


2.636( 


1.9 


K\ 


: sU 2 3 P 4 


2.490( 20.0) 


2.757( 


1.6 


0: 


ss l 3 Si 


1.019( 0.0) 


1.033( 


2.2 




ss 1 3 P 


1.370( 40.0) 


1.319( 


0.0 


h ■ 


SS l 3 Pi 


1.512( 4.0) 


1.533( 


0.3 




ss 1 3 P 2 


1.525( 5.0) 


1.493( 


0.3 


0: 


ss 2 3 Si 


1.680( 20.0) 


1.850( 


0.8 


0: 


ss 1 3 P 3 


1.854( 7.0) 


1.848( 


0.0 


/ 2 : 


ss 2 3 P 2 


2.011( 69.0) 


2.160( 


0.1 


h- 


ss 3 3 P 2 


2.297( 28.0) 


2.629( 


1.6 


n : 


ud l 1 ^ 


0.140( 0.0) 


0.144( 


0.2 


P ■ 




0.767( 1.2) 


0.792( 


0.1 



16 



h 


: ud l l P x 


1.231( 10.0) 


1.392( 


2.1) 


a 


: ud 1 3 P 


1.450( 40.0) 


1.491( 


0.0) 


ai 


: ud l 3 Pi 


1.230( 40.0) 


1.568( 


0.7) 


a 2 


: ud 1 3 P 2 


1.318( 0.7) 


1.310( 


0.0) 


7r : 


2 1 S'o 


1.300( 100.0) 


1.536( 


0.1) 


P ■ 


ud 2 3 5i 


1.465( 25.0) 


1.775( 


1.4) 




: l 1 /^ 


1.670( 20.0) 


1.870( 


0.9) 


P ■ 


ud l 3 £>i 


1.700( 20.0) 


1.986( 


1.9) 


P3 


: wrf 1 3 D 3 


1.691( 5.0) 


1.710( 


0.0) 


7r : 


Mrf 3^0 


1.795( 10.0) 


2.166( 


7.9) 


p ■ 


m3 S 3 ^ 


2.149( 17.0) 


2.333( 


0.7) 


pi 


: ud 1 3 F 4 


2.037( 26.0) 


2.033( 


0.0) 


7T 2 


: m3 2 1 D 2 


2.090( 29.0) 


2.367( 


0.5) 


P3 


: ud 2 3 D 3 


2.250( 45.0) 


2.305( 


0.0) 


P5 


: 1 3 G 5 


2.330( 35.0) 


2.307( 


0.0) 


P6 


: ud l 3 H 6 


2.450( 130.0) 


2.547( 


0.0) 



101.0 



(The numbers in parentheses represent experimental uncertainties and x 2 contributions 
for each meson.) 

The heavy upsilon fits to the ground and excited states are due more to the specific 
form of the Adler-Piran potential than to the relativistic features of the equations, since 
the heavy quark motions are nearly nonrelativistic. The spin-orbit splittings arise from our 
apportionment of the Adler potential into scalar and vector parts. The good fits to the heavy 
B mesons are due primarily to the fact that in the infinite quark mass limit our equations 
reduce to the one body Dirac equation. 

For charmonium, the semirelativistic effects of the formalism become important, the 
higher order relativistic effects becoming more important in the D and then mesons. The 
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only significant weakness in our fits (due to limitations of apportioning our potential only 
between scalar and vector potentials) manifests itself at the LS multiplet level for the lighter 
mesons. For those states, partial inversion takes place. 

However, the equations work well in the extreme relativistic domain of the very light 
mesons and their lower excitations - (K,K*,n,p) and first radial excitations. 

In our equations, the pion is a Goldstone boson in the sense that its mass tends toward 
zero numerically in the limit in which the quark mass numerically goes toward zero. This 
may be seen in the accompanying plot. Note that the p meson mass approaches a finite 
value in the chiral limit. This also holds true for the excited pion states. 

pi and rho masses vs. quark mass 

800 
700 
600 
500 
400 
300 
200 
100 


5 10 15 20 25 30 35 40 45 50 55 

Fig. 1.7T and p masses versus quark mass (in MeV) 
Our results have room for improvement in the spin-orbit splittings of the light mesons 
and other orbital and radial excitations. This is a consequence of our assumption that at 
long distance the confining interaction is pure scalar. 

With just two parametric functions, (amounting to an ad hoc division of the Adler 
potential into vector and scalar parts), we are able to obtain a fit about as good as that 
obtained by Godfrey and Isgur |2T|, who use six parametric functions, basically one for each 
type of spin dependence. The Godfrey and Isgur results seem more accurate than ours for 
the spin-orbit splitting of the light mesons, likely due to their use of different parameters for 
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each spin-dependent part. Note that they must also introduce parameters that smooth the 
potential, a procedure that is not necessary in the two-body Dirac equations because of the 
connection of their structure to that of the one-body Dirac equation. 



IV. TWO-BODY DIRAC EQUATIONS FOR GENERAL CO VARIANT 

INTERACTIONS 

One may generalize the interactions in the Two-Body Dirac Equations to ones other 
than vector and scalar. How does one determine the forms of these equations for arbitrary 
interactions? We use this opportunity to discuss the role played by supersymmetry in the 
origins of the constraint formalism for two-body Dirac equations [|K],|ll[ • 



A. The Role of Supersymmetry 

We first rewrite one-body Dirac equation in terms of theta matrices defined by 

^ = ^757^ = 0,1,2,3, # 5 = ^75, (4.0a) 

r,n + = -?r, [0 5 ,0"] + = o, [0 5 ,0 5 ]+ = -i. (4.0b) 

In terms of these the free Dirac equation and its Klein-Gordon square become 

Sip = (p- 9 + m9 5 )4i = (4.1) 

in which S is literally the "operator square root" of Einstein's mass shell condition: 

S 2 ip = ~(p 2 + m 2 )ip = 0. (4.2) 

The supersymmetry is most easily revealed through examination of "pseudoclassical" me- 
chanics - the "correspondence-principle" limit of the Dirac equation. In that limit, these 
theta matrices become Grassmann variables obeying "pseudoclassical" Berezin-Marinov 
brackets which upon quantization generate the Dirac algebra |21|: 6 U } = zj/*", {#5, 6^} = 

o,{0 5 ,0 5 } = i. 
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In this limit, the Dirac equation becomes a constraint imposed on both bosonic (p) and 
fermionic (9, 9 5 ) variables: 

S = (p ■ 6 + m6 5 ) « (4.3) 

while the mass shell condition comes from the bracket 

H = -{S,S} =p 2 + m 2 « 0. (4.4) 

We first examine the supersymmetry of the free Dirac constraint under transformations 
generated by the self-abelian Q defined by 



g = p ■ 6 + ^f-p 2 6 5 , (4.5a) 

{G,G} = 0. (4.5b) 

That is, 

{G,S}^0. (4.6) 

Q is a supersymmetry generator in the sense that it mixes bosonic variables like the mo- 
mentum with the Grassmann or fermionic variables. To maintain this supersymmetry in 
the presence of interaction, we must introduce interactions that have zero brackets with 
Q. The position four-vector x is not supersymmetric and in fact displays pseudoclassical 
zitterbewegegung. However, the "zitterbewegungless" position variable 

x» = ^ + ^ (4.7) 
m 

is supersymmetric. But, because of the presence of m in this variable, this object must 
itself be modified in the presence of scalar interaction M = m + S. We accomplish this by 
generalizing this position variable to the self-referent form 

x» = x » + JL*. (4.8) 
M(x) v 1 

20 



This is a supersymmetric position variable appropriate for scalar interactions in the sense 
that 

{Q,x}^0. (4.9) 
The supersymmetric constraints (both fermionic and bosonic) then become 

S = p ■ 9 + M(x)9 5 « 0, -{S, S} = H = p 2 + M 2 (x) « 0. (4.10) 
Since #5 = ,the expansion of the self-referent form truncates 

M(x)=M(x)+ (4.11) 
so that the S constraint assumes the standard form 

S = p-6 + M(x)6 5 ^0 (4.12) 
while the mass-shell constraint becomes the familiar square 

-{s, s} = n = p 2 + m 2 {x) = P 2 + m 2 {x) + 2idM(x) ■ ee 5 . (4.13) 

i 

One arrives back at the Dirac equation and its square for scalar interaction by replacing 
Grassmann variables with theta matrices and dynamical variables x and p with their operator 
forms. 

SiIj = [p ■ 6 + M(x)6 5 ]iIj = (4.14a) 

Hip = [p 2 + M 2 (x) + 2idM(x) ■ ee 5 ]ip = 0. (4.14b) 

These contain the usual spin-dependent corrections expected for scalar interactions. We 
have gone through this exercise in order to show that this structure is present in the usual 
result. The supersymmetry generated by Q and realized through the presence of a; is a 
natural feature of both the free Dirac equation and its standard form for external scalar 
interaction. 
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How does one implement this supersymmetry in the case of two interacting particles? 
For two pseudoclassical free particles we begin with 



Sio = Pi ■ 9i + mi9 5i w 0, % = 1, 2 (4.15) 

in which the two sets of 0's are independent Grassmann variables so that the mutual S i0 
bracket vanishes strongly: 

{S w ,S 20 } = 0. (4.16) 

In the presence of interaction, we require the preservation of supersymmetry for each spinning 
particle. For scalar interactions we accomplish this through 

m t -> Mi( Xl - x 2 ) -> Mi(xi - x 2 ) = Mi, i = l,2 (4.17) 

depending on a supersymmetric variable for each particle position 

# = tf + ^i, < = 1,2. (4.18) 

Note that the Grassmann Taylor expansions of the Mj truncate. Carrying out those expan- 
sions, we obtain the following two-body Dirac constraints: 

d M 9 9 9 

S 1 = Pl ■ 9 1 + M,9 5l =p 1 -6 1 + M 1 9 51 - i 1 \l 52 51 » 0, (4.19a) 

S 2 =p 2 -9 2 + M 2 9 52 =p 2 -9 2 + M 2 9 52 + % 2 - 2 51 52 w 0. (4.19b) 

whose "squares" become 

-^S^^H^vl + Ml 2 = 1,2. (4.20) 

One can show that all four S and 7Y's are mutually compatible provided that two other 
conditions are met. Straightforward computation leads to 

{5i,5 2 } = -Pi • ^9 51 9 52 -p 2 ■ ^9 51 9 52 . (4.21) 

22 



This bracket vanishes strongly provided that 

d{Ml - Ml) = 0, (4.22) 
(the relativistic "third law" condition), and 

M { = Mi(x±) (4.23) 

in which 

x ± = (pi +p 2 ) 2 ~ 

which, eliminates covariantly the troublesome "relative time". Now we use these results to 
rewrite the constraints in a form which we can easily generalize to interactions other than 
scalar. The third law condition Eq. ( |4.22j ) has the solution 



Ml - Ml = m\-m\ (4.25) 
which has the convenient invariant hyperbolic parametrization 

Mi = rrii chL + m 2 shL, M 2 = m 2 chL + mi shL, (4.26) 

in which 

L = L(x ± ). (4.27) 

In this case just one invariant function suffices to define the two-body interaction. For such 
interactions 

U x - H 2 = v\ + M\ -p\ - Ml = p\ + m{ -p\ - ml w (4.28) 



just as in the noninteracting case (and indeed also in the spinless interacting case [p0|,|T0|] ) . 
When the potentials satisfy these conditions, the Dirac-like constraints S\ and S 2 are com- 
patible and generate a consistent "pseudoclassical" dynamics. Note that the combination of 
the constraints corresponding to the difference of the squares eliminates the relative energy 
in the cm. rest frame. To see this write the constraints in terms of the total four-momentum 
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pM= p M +p M. _p2 _ w 2. p» = pM/ W; (429) 

and relative four-momentum 

P' 1 = (^2 - t\V2)l w \ e i + e 2 = «>, ei - e 2 = {m\ - m 2 2 )/w. (4.30) 

(Note that p is canonically conjugate to x±. The e^'s are the invariant cm. energies of each 
of the (interacting) particles.) Thus the constraint difference 

Hx -H 2 = -2P-p « (4.31) 

places an invariant restriction on the relative energy. The other linearly independent com- 
bination 

(e 2 Wi + e{H 2 )/w = p\ + <& w - b 2 {w) « 0, (4.32) 

in which 

6 2 ( w ) = (u, 4 - 2w 2 (m? + m 2 2 ) + {ml - m 2 2 ) 2 )/4w 2 (4.33) 

incorporates exact relativistic two-body kinematics and governs the dynamics through the 
quasipotential 

$ w = M 2 -m\ = Ml - m\. (4.34) 

Just as in the one-body case, we canonically quantize this system by replacing the Grassmann 
variables 9^, 9 5i i — 1, 2 by two mutually commuting sets of theta matrices, and 

{x",p v } -> [x",p v ] =irf v . (4.35) 

This leads to strongly compatible ([«Si,«S 2 ] = 0) two-body Dirac equations (for scalar inter- 
actions) in minimal interaction form 

Sii) = (0i • p + eA • P + M,e 51 - idL ■ 9 2 9 52 9 51 )ijj = 0, (4.36a) 

s 2 ^j = {-e 2 ■ p + e 2 e 2 ■ P + M 2 e 52 + idL ■ e l e b2 e bl )^ = o, (4.36b) 
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in which 



dL = 



dM l 



dM 2 



(4.37) 



M 2 



M 1 ' 



Mi = mi chL + m 2 shL, M 2 = m 2 chL + mi shL. 



(4.38) 



We see that in these coupled Dirac equations the remnants of pseudoclassical supersymme- 
tries are the extra spin dependent recoil corrections to the ordinary one-body Dirac equa- 
tions. Without those terms the two equations would not be compatible. When Mj = rrii + Si, 
these extra terms vanish when one of the particles becomes infinitely heavy. 

B. Hyperbolic form of the two-body Dirac equations for General Covariant 

Interactions 

How do we introduce general interactions? We accomplish this by recasting the mini- 
mal interaction forms of the two-body Dirac equations into a more general form, one that 
generalizes the hyperbolic forms we encountered above. Simple identities such as 



ch 2 (A) - sh 2 (A) = 1 



(4.39) 



lead to 



= (ch(A)S 1 + s/i(A)S 2 )V = 0, 



(4.40a) 




(4.40b) 



in which appear auxiliary constraints defined by 



= (S 10 ch(A) + S 20 sh(A))i> = 0, 



(4.41a) 



with 



S 2 ^ = (S 20 ch(A) + S w sh(A))rl> = 0, 
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(4.41b) 



A = -e hl e h2 L{ Xl _). (4.42) 

Note that in this form, the interaction enters only through an invariant matrix function A 
with all other spin-dependence contained in the kinetic free Dirac operators <Sio,c>2o- One 
can show that both Si and Sj constraints are compatible for general A: 



[Si, S 2 ]^ = and [S u S 2 ]i> = (4.43) 

provided only that 

A = A(zj_). (4.44) 

For the polar interactions we find 

A(xj_) = —L(xj_)9 5 i9 52 scalar (4.45a) 

A(x ± ) = J(x±)P ■ 9 X P ■ 9 2 time like vector (4.45b) 

A(xj_) = Q{x±)9i± ■ 9 21 _ space like vector (4.45c) 

A(arj.) = F(x±)0 1± ■ 9 2± 9 51 9 52 P ■ 9 X P ■ 9 2 tensor (polar). (4.45d) 

The constraint equations for QED presented at the beginning of this paper are generated in 
this form by taking L = = T and Q = —J = — \ln(l — 2A/w). For the QCD calculations 
which we performed above L is non zero and determined [|32| by the S of Eq. (|3.5b| ) while 
the vector interaction is written in terms of A of Eq. (|3.5a|) and just as for QED has the 



Feynman gauge combination J = —Q. For the axial counterparts the hyperbolic forms of 
our constraints are (note the minus sign) 

Slip = (cfc(A)Si - sh(A)S 2 )tfj = (4.46a) 



Slip = (ch(A)S 2 - s/i(A)Si)V> = 0, (4.46b) 
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in which Si and S2 are defined as before while the interactions appear through 

A(x±) = C(x±)/2 pseudoscalar 



(4.47a) 



A(x±) = H(x±)P ■ Q\P • 2 5 i9 52 time like pseudovector 



(4.47b) 



A(x±) = I(x±)9i± ■ 9 2 ±0 51 9 52 space like pseudovector 



(4.47c) 



A(x±) = Y(x±)9 1± ■ 9 2± P ■ QiP ■ 9 2 tensor (axial). 



(4.47d) 



Future research will determine the relative importance of the interactions other than scalar 
and vector in meson spectroscopy. 

We conclude by listing important features realized by these two body Dirac equations 
which are incompletely realized if present at all in other approaches to the relativistic two- 
body problem. A) manifest covariance, B) simple three-dimensional Schrodinger-like forms 
similar to their nonrelativistic counterparts, C) spin dependence dictated by the Dirac-like 
structure of the equations (not put in by hand), D) thoroughly tested in QED: the equations 
reproduce the correct perturbative results when solved nonperturbatively (they qualify as 
bona fide wave equations less likely to produce spurious nonpeturbative effects), E) the close 
connection to the one body Dirac equations automatically eliminates the ad hoc introduction 
of cutoff parameters, G) its potentials are determined through straightforward connection 
to perturbative quantum field theory or introduced semiphenomenologically. 
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